Photon storage in A-type optically dense atomic media. III. Effects of inhomogeneous 

broadening 



O 

o 

(N 

a; 

00 
00 

(N 
> 

oo 
o 

(N 
O 

9 L,: 
^— > ■ 

c ■ 

s : 
cr 



X 



Alexey V. Gorshkov, 1 Axel Andre, 1 Mikhail D. Lukin, 1 and Anders S. S0rensen 2 

1 Physics Department, Harvard University, Cambridge, Massachusetts 02138, USA 
2 QUANT OP , Danish National Research Foundation Centre of Quantum Optics, 
Niels Bohr Institute, DK-2100 Copenhagen 0, Denmark 
(Dated: February 1, 2008) 

In a recent paper [Gorshkov et ai, Phys. Rev. Lett. 98, 123601 (2007)] and in the two preceding 
papers [Gorshkov et ai, Phys. Rev. A 76, 033804 (2007); 76, 033805 (2007)], we used a universal 
physical picture to optimize and demonstrate equivalence between a wide range of techniques for 
storage and retrieval of photon wave packets in homogeneously broadened A-type atomic media, 
including the adiabatic reduction of the photon group velocity, pulse-propagation control via off- 
resonant Raman techniques, and photon-echo-based techniques. In the present paper, we generalize 
this treatment to include inhomogeneous broadening. In particular, we consider the case of Doppler- 
broadened atoms and assume that there is a negligible difference between the Doppler shifts of the 
two optical transitions. In this situation, we show that, at high enough optical depth, all atoms 
contribute coherently to the storage process as if the medium were homogeneously broadened. We 
also discuss the effects of inhomogeneous broadening in solid state samples. In this context, we 
discuss the advantages and limitations of reversing the inhomogeneous broadening during the storage 
time, as well as suggest a way for achieving high efficiencies with a nonreversible inhomogeneous 
profile. 

PACS numbers: 42.50.Gy, 03.67.-a, 32.80.Qk, 42.50.Fx 



I. INTRODUCTION 

The faithful storage of a traveling light pulse in an 
atomic memory and the subsequent retrieval of the state 
are currently being pursued in a number of laboratories 
around the world. A central question that emerges from 
these studies is which approach represents the best pos- 
sible strategy and how the maximum efficiency can be 
achieved. In a recent paper we used a novel univer- 
sal physical picture to optimize and demonstrate equiv- 
alence between a wide range of techniques for storage 
and retrieval of photon wave packets in A-type atomic 
media in free space, including the adiabatic reduction of 
the photon group velocity, pulse-propagation control via 
off-resonant Raman techniques, and photon-echo-based 
techniques. In two preceding papers Refs. [HQ j which 
we will refer to henceforth as papers I and II, we present 
the full details of the analysis of Ref. [l[ as well as sev- 
eral of its extensions, both for an ensemble enclosed in a 
cavity and for the free-space model. While the analysis 
of papers I and II is limited to homogeneously broadened 
media, many experimental realizations, such as the ones 
using warm atomic vapors [J1 or the ones using impu- 
rities in solid state samples [a, Q , have some degree of 
inhomogeneous broadening, whose presence will modify 
the optimal control strategy and the values for the max- 
imum efficiency. The subject of the present paper is the 
extension of the analysis of papers I and II to include 
inhomogeneous broadening. 

The remainder of the present paper is organized as fol- 
lows. In Sec. HH we discuss the effects of inhomogeneous 
broadening assuming that the atoms fully redistribute 
themselves between frequency classes during the storage 



time, which would be the case, for example, in Doppler- 
broadened atoms with sufficiently long storage times. In 
particular, we optimize the storage process and show that 
at high enough optical depth, all atoms contribute coher- 
ently as if the medium were homogeneously broadened. 
Then in Sec. IIII1 we discuss the effects of inhomogeneous 
broadening assuming there is no redistribution between 
frequency classes during the storage time, which would 
be the case in atomic vapors for short storage times or 
in solid state samples. In particular, we discuss the ad- 
vantages and limitations of reversing the inhomogeneous 
broadening during the storage time [7| , as well as suggest 
a method for achieving high efficiencies with a nonre- 
versible spectrally well-localized inhomogeneous profile. 
In Sec. IIV[ we summarize our analysis of the effects of 
inhomogeneous broadening. Finally, in the Appendix, we 
present some details omitted in the main text. 



II. INHOMOGENEOUS BROADENING WITH 
REDISTRIBUTION BETWEEN FREQUENCY 
CLASSES DURING THE STORAGE TIME 



In this section, we consider a particular case of inho- 
mogeneously broadened media: the case of a Doppler- 
broadened atomic vapor in free space. We first describe 
our model in Sec. Ill AI We then use this model in 
Sec. Ill Bl to analyze storage and retrieval of photons in 
Doppler-broadened media. 
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FIG. 1: (Color online) A-type medium coupled to a classical 
field (solid) with Rabi frequency Q(t) and a quantum field 
(dashed) . Collective enhancement |J effectively increases the 
atom-field coupling constant g up to gvN, where N is the 
number of atoms in the medium. Aj is a frequency shift due 
to inhomogeneous broadening. 



A. Model 

As in paper II, we consider a free-space medium of 
length L and cross-section area A containing N — 
J dzn(z) atoms, where n(z) is the number of atoms 
per unit length. We assume that within the interac- 
tion volume the concentration of atoms is uniform in the 
transverse direction. The atoms have the A-type level 
configuration shown in Fig. [TJ They are coupled with a 
collectively enhanced coupling constant gvN (g assumed 
to be real for simplicity) to a narrowband quantum field 
centered at a frequency lo\ and described by a slowly 
varying operator £{z,t). The atoms are also coupled to 
a copropagating narrowband classical control field cen- 
tered at frequency L02 with a Rabi frequency envelope 
t) = il(t — z/c). We assume that quantum electro- 
magnetic field modes with a single transverse profile are 
excited. As discussed in detail in paper II, we neglect 
reabsorption of spontaneously emitted photons and treat 
the problem in a one-dimensional approximation. 

In order to model the inhomogeneous broadening, we 
sort all the atoms into separate velocity classes labeled 
by j, such that all atoms in a certain velocity class have 
approximately the same velocity Vj . We will assume that 
during the experiment the atoms do not change their po- 
sitions significantly compared to the longitudinal varia- 
tion of the fields. Therefore, the Doppler effect will be 
the only effect of nonzero atomic velocities that we will 
consider. 

Furthermore, we will assume that the difference in the 
Doppler shifts of the two optical transitions can be ne- 
glected so that atoms in all velocity classes can simulta- 
neously stay in two-photon resonance with the two fields. 
For copropagating beams, this assumption is exactly sat- 
isfied if the two metastable states are degenerate. More- 
over, below we shall consider a storage technique which 
we refer to as fast storage [l|, H| , where the control field 
consists of a simple 7r pulse, which will work perfectly pro- 



vided its Rabi frequency is much greater than the broad- 
ened linewidth. In this case, the assumption of equal 
Doppler shifts of the two optical transitions is thus not 
needed. However, we shall also consider the so-called adi- 
abatic storage schemes 0, Q , where the splitting of the 
two metastable levels by a nonzero frequency difference 
uj sg , as well as Doppler broadeningoccurring for nonco- 
propagating fields, will play a role [8|. In order to ensure 
that the difference in Doppler shifts has a negligible ef- 
fect in this situation, we consider copropagating beams 
and assume that the total accumulated phase difference 
Tvuj sg /c is much less than unity, where v is the ther- 
mal speed of the atoms, c is the speed of light, and T is 
the duration of the incoming quantum light pulse. This 
condition is usually satisfied even in room-temperature 
experiments, such as the experiment using Rb vapor 
in Ref. Q, where uj sg = (2tt)6.8 GHz and T ~ 200 fjs, 
which gives Tvuj sg /c < 0.01. 

For simplicity, we also assume that the velocities do not 
change during the processes of storage and retrieval but 
fully rethermalize during the storage time [T, T T ] . The 
atoms thus fully redistribute themselves among different 
velocity classes during the storage time, and at time T x 
the spin wave is the same across all velocity classes. In 
addition to being relevant for Doppler broadening, much 
of the discussion of this section will also apply to solid 
state systems with inhomogeneous broadening. The fact 
that solid state impurities do not have the redistribution 
between frequency classes during the storage time (which 
we assume in this section) will, however, introduce mod- 
ifications. Some of these modification are discussed in 

sec. inn 

We assume that the control is detuned by A with 
respect to stationary atoms, 0J2 = oj es — A, while the 
quantum field is in two-photon resonance, i.e., u)\ = 
Lu eg — A (where oj es and uj eg are atomic transition fre- 
quencies). We define the same slowly varying operators 
as in Eqs. (A7)-(A11) in paper II, except now, at each 
z and t, we have continuous atomic operators for each 
velocity class. For example, we define 

N i 

ivz '(=1 

where j indicates the velocity class, i runs over the N° z 
atoms near z in the velocity class j, pj is the fraction 

of atoms in the velocity class j, YljPj = 1, and ajiu = 
\(J>)jU) {v\ indicates the internal state operator between 
states and \v) for the ith atom in the jth velocity 
class. 

As in papers I and II, we use the dipole and rotating- 
wave approximations, assume that almost all atoms are 
in the ground state at all times, and consider equations 
of motion only to first order in £. As described in pa- 
pers I and II, under reasonable experimental conditions, 
the incoming noise corresponding to the decay of atomic 
coherences is vacuum, and efficiency is the only number 
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we need in order to fully characterize the mapping. Fol- 
lowing then the same steps as in Sec. II of paper II, we 
obtain complex number equations 



(fit + cd z )£ = igVNj2 VFiPjn[z)L/N, 



dtP, 

dtS 3 = iQ*Pj 



(2) 



- (j+i ( A+ Aj) )Pj +igVN^/p]£ +i£lSj,(3) 

(4) 



■U) e 



where 7 is the polarization decay rate, Aj 
to es Vj I c is the Doppler shift due to the velocity Vj of the 
jth velocity class, Pj is the complex number representing 
the optical polarization \fNa ge , and Sj is the complex 

number representing the spin wave VNa gs . We assume 
in Eq. (TJ| that the decay rate of the spin wave Sj is neg- 
ligible. As in Eqs. (5)- (7) in paper II, we now go into the 
comoving frame t' — t — z/c, introduce the dimension- 
less time t = jt' and a dimensionless res caled co ordinate 
z = J Q dz'n(z')/N , absorb a factor of yjc/jlsy) into the 
definition of £, and obtain 



(5) 



d- t Pj = -(1 + i(A + Aj))Pj + iVd^p-£ + iQSj, (6) 

d- t Sj = m*Pj, (7) 

where we have identified the optical depth d = 
g 2 NL/(jc), and where Q = £1/7, A — A/7, and 
Aj = Aj/7. Note that d is defined here as the opti- 
cal depth that the sample would have had if there were 
the same number of atoms but no inhomogeneous broad- 
ening. This quantity should not be confused with the ac- 
tually measured optical depth d! . With inhomogeneous 
broadening, the measured value will be lower: d' < d. 
Later, we shall derive explicit relations between these two 
quantities. It is essential to realize that both quantities 
d and d' play a role in the performance of the memory, 
as we shall see below. 

We assume that all atoms are initially pumped into the 
ground state, i.e., no P or S excitations are present in 
the atoms. We also assume that the only input field ex- 
citations initially present are in the quantum field mode 
with a normalized envelope shape £- m (t) nonzero on [0, T] 
(where T = Tj). The goal is to store the state of this 
mode into some spin-wave mode and at a time T r > T 
retrieve it onto an output field mode. The initial condi- 
tions for storage are Sj(z, 0) = and Pj(z, 0) = for all 
j and £(0,t) = £in(t), and the storage efficiency is given 
by 

(number of stored excitations) f 1 ^,,2 

Vs = ~, r 7^ : 7— r=/ dz\S{z,T)\ , 

(number ot incoming photons) Jq 

(8) 

where S(z,T) — y/p]Sj ;(z, T) is the spin wave, to 
which all Sj average after rethermalization. The ini- 
tial conditions for retrieval are £(0,t) — 0, and, for 



Sj(z, T r ) = y/pjS(l — z,T) for forward or backward re- 
trieval, respectively (as in paper II). The total efficiency 
of storage followed by retrieval is then given by 



Vtot 



(number of retrieved photons) 
(number of incoming photons) 



dt\£ ont (t)f 



(9) 

where £ ou t(t) = £(l,t). If during retrieval we instead 
start with Sj(z,T T ) = ^/pjS(z) for some normalized spin 
wave S(z), then the same equation will give the retrieval 
efficiency from this mode: 



Vr 



(number of retrieved photons) 
(number of stored excitations) 



dt\£ out (t)\< 



(10) 

For completeness, we note that, for the cavity model 
described in paper I, the equations corresponding to 
Eqs. ([5j)- (JTJ) above are (without rescaling of variables) 



^out 

dtPj 
d t Sj 



£ in + i^Cj2VPj p j> (11) 

j 

-(7 + i(A + Aj))Pj - iC^/pjJ2 VP^ p k 



+iQ,Sj + iy / 2-fC^/p-£ in , 



ifL*Pi 



(12) 
(13) 



all j, Pj(z,T r ) = 0, and Sj{z,T T ] 



/p]S(z,T) or 



where C is the cooperativity parameter equal to the opti- 
cal depth of the ensemble multiplied by the cavity finesse. 
The initial conditions for storage in the cavity model are 
Sj(0) = and Pj(0) = for all j and £ in (t) ^ 0, while 
the initial conditions for retrieval are Sj(T r ) = ^/pjS(T) 
and Pj(T T ) = for all j and £i n {t) = 0. It is assumed 
that during the storage time atoms rethermalize and all 
Sj average to the same value S(T) = J2j ^/PjSjiT). 

In the homogeneously broadened case discussed in pa- 
pers I and II, we defined the so-called adiabatic and fast 
regimes for storage and retrieval. Both of these limits 
can also be achieved in the presence of Doppler broaden- 
ing. The adiabatic regime corresponds to smooth control 
and input fields such that the optical polarization Pj in 
Eq. (J6j> (or Eq. (fT2|) for the case of the cavity model) can 
be adiabatically eliminated. The fast regime corresponds 
to storage and retrieval techniques in which the control 
field consists of a very short and powerful resonant it 
pulse between states |s) and |e). The only difference 
in the requirements from the homogeneously broadened 
case is that now the control field SI in the tt pulse must 
also be much greater than the inhomogeneous width. We 
refer the reader to paper I for a detailed discussion of the 
adiabatic and fast photon storage techniques, as well as 
for a full list of references. 

Equations ©-([7]) can be solved numerically by intro- 
ducing sufficiently many discrete velocity classes. How- 
ever, when the control is constant (i.e., a step) or when we 
are in the fast limit, these equations can be solved with- 
out discretizing the velocity distribution by using Laplace 
transformation in time, t — > v. In this case, the inverse 
Laplace transform has to be taken numerically at the end. 
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Alternatively, in the case of retrieval alone or in the case 
of storage followed by retrieval, if one is interested only 
in the efficiency and not in the output mode, one can 
compute this efficiency both in the free-space model and 
in the cavity model without computing £ ou t(t) directly 
from the Laplace transform £ on t(v) as 



1 r°° 

Vr = 2^ J d£\£ out (v = 



(14) 



Below we shall use both the numerical method with dis- 
crete velocity classes and the method of the Laplace 
transformation in time. 



An important result in the discussion of homoge- 
neously broadened ensembles in papers I and II was the 
result that the retrieval efficiency is independent of the 
shape of the control and the detuning, provided that all 
excitations are pumped out of the system. Moreover, in 
both the cavity and the free-space cases we were able 
to deduce in papers I and II explicit formulas for the 
retrieval efficiency. Although in the inhomogeneously 
broadened case discussed in the present paper we have 
not been able to find an explicit formula for the retrieval 
efficiency, we will present now the proof that even with 
inhomogeneous broadening the retrieval efficiency is in- 
dependent of the detuning A and the control shape fi(£). 

We consider first the cavity model given in Eqs. ([TT]) - 
(I13[) . Since we are interested in retrieval, we set £ ln = 0. 
We also consider a general situation, in which Sj(t = 0) 
are not necessarily equal for different velocity classes j 
(we also shifted for simplicity the beginning of retrieval 
from t = T T to t = 0). Then, using Eq. (fTTj) . the retrieval 
efficiency is 



/>oo t>OC 

rk = / dt \£ out (t)\ 2 = 2 7 C]T ^pjp~k / dtP 3 {t)P* k {t), 
Jo j>k Jo 

(15) 

and the following identity can be explicitly verified from 
Eqs. CUD and fB} : 



S-j SV , 



-(2 7 + i(A j -A k ))P j P k * (16) 



If M is the number of velocity classes, Eq. ([16]) stands for 
M 2 equations in M 2 variables PiP k - We can write them 
in matrix form and, in principle, invert the M 2 x M 
matrix on the right-hand side and thus solve for PjP k as 
a linear combination of ^{PaP^ + S a S£) for various a 
and b. Inserting this into Eq. (fT5|) . applying the funda- 
mental theorem of calculus, and assuming the retrieval is 
complete (i.e., no excitations remain in the atoms), the 
retrieval efficiency can be expressed as a linear combi- 
nation of S a (0)«S^(0), and is thus independent of control 
and detuning. 



In Appendix we present an analogous derivation, 
which shows that the free-space retrieval efficiency in the 
presence of inhomogeneous broadening is also indepen- 
dent of detuning and control. Numerical calculations also 
show that adiabatic elimination of Pj , as in the homoge- 
neously broadened case discussed in papers I and II, does 
not change the exact value of the efficiency. 

Since the retrieval efficiency is thus independent of the 
exact method used for retrieval, we shall here mainly con- 
sider the fast retrieval from S(z). In this case Laplace 
transformation in time can be used to solve the problem 
analytically. We will focus for the rest of this section 
on the free-space model. We assume that the retrieval 
tt pulse arrives at t = and that it perfectly trans- 



B. Retrieval and storage with Doppler broadening fers Pj(z) 



S,(z) 



0. 



0, Sj(z) = vpj5(z) to Pjii) = i^pjS(z), 
After the tt pulse, we Laplace-transform 



Eqs. ([5]) and ([6]) in time t — > v and obtain 



dz£ = iVd^y^^ 



PJP; 



(17) 



VP, 



/pJS(z) = -(l + iA^Pj+iVd^/pJS. (18) 

Solving for Pj from the second equation and inserting it 
into the first equation, we find 



£{z = l,v) 



dzS{S)Vdf(v)e- df{v){1 - l \ (19) 



where 



1 



+ v + iAj 



dAp(A)- 



1 



1 + v + iA 

(20) 

and where p(A) is the Doppler profile. 

In Doppler-broadened media, the resonant optical 
depth is reduced by a factor of ~ j/Aj, where Ai is 
the width of the (inhomogeneous) Doppler profile. The 
naive expectation would therefore be that we could sim- 
ply treat Doppler-broadened atoms as Doppler-free but 
with a reduced optical depth. As we will show below, this 
prescription would be correct if the broadened line shape 
were a Lorentzian, as considered, e.g., in Ref. [io| . For a 
Gaussian profile this prescription is, however, not appli- 
cable, and qualitatively different behavior is obtained. 

To proceed, we first evaluate f(v) for three different in- 
homogeneous profiles p(A). For a homogeneously broad- 
ened ensemble, the line shape function p(A) is just a S 
function: 



p(A) = S(A), 

m 1 



1 + v 



(21) 
(22) 



If we have a Lorentzian inhomogeneous profile with Aj 
half width at half maximum (HWHM), we get (with Ai = 
A1/7) 



p(A) 



A, 



1 



/(«) = 



TT A 2 + A 2 

1 



l + v + Ai 



(23) 
(24) 
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For a Gaussian inhomogeneous profile with (rescaled by 
7) standard deviation a (and rescaled HWHM Ai = 
oV21n2), 



P (A) 

M = 



V2~7 



: exp 



2 CT 2 



e 2 ° 




(25) 
(26) 



where the last equality assumes Re[v) > 0, and the 
complementary error function is defined as erfc(x) = 
1 — 7r _1//2 2 Jq exp(— x' 2 )dx' . Using the definition of erfc, 
we can analytically continue f(v) into Re[u] < 0. 

If we insert the Lorentzian result for f(v) from Eq. 
into Eq. (fl"9]l and rescale v, we find that, compared to ho- 
mogeneous broadening (Eq. ([22]) ). the Lorentzian broad- 
ening effectively just replaces 7 with 7 + Ai, which is 
equivalent to reducing d by a factor of 1/(1 + A1/7) 
(which for Ai ^> 7 is equal to 7/A1). Therefore, since we 
have shown that the retrieval efficiency is independent of 
the retrieval method, the naive rescaling of d to dj/Aj 
can indeed be used to calculate the retrieval efficiency 
when the broadening is Lorentzian. Similarly, the same 
can be shown for the cavity model. 

For Doppler broadening, p(A) is, however, Gaussian 
as in Eq. (|25|) . Using Eqs. (fT4|) and (fl9l) . we can write 
the retrieval efficiency in the form 



r,[S(z)] 



dz 



dz k(z, z )S(z)S*(z ) 



(27) 



for some complicated kernel k. Applying the iterative 
technique used in Sec. Ill of paper II or by directly diago- 
nalizing k on a grid, we can compute the optimal retrieval 
modes (i.e., the eigenvectors with the largest eigenvalues) 
for each d and a. 

Before plotting and analyzing the optimal spin waves 
and the maximum efficiency obtained using Eq. (|27[) . let 
us discuss what we expect. Assuming Ai 3> 7, the reso- 
nant optical depth is reduced to d! — d(-f / Ai)v V ln2 in 
the presence of Gaussian broadening. In contrast to re- 
trieval with a Lorentzian profile, however, retrieval with 
a Gaussian profile is not equivalent to Doppler-free re- 
trieval with reduced optical depth (as we can see by com- 
paring Eqs. (|26p and (|22|) L Moreover, we will show now 
that, with true (Gaussian) Doppler broadening, at high 
enough optical depth all atoms contribute coherently as 
if the medium were homogeneously broadened, which is 
the main result of this section. Although this result holds 
for any control, it is most easily explained in the case of 
fast retrieval: after the tt pulse, the spontaneous emission 
(or more precisely the polarization decay at a rate 7) and 
the dephasing due to the inhomogeneous broadening will 
cause the polarization P(t) = • y/pjPj(t) (with an ini- 
tial velocity-symmetric polarization Pj(0) = y/pJP(0)) 
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FIG. 2: Optimal spin- wave modes to retrieve from (in the 
forward direction) at the indicated values of d' in the presence 
of (Gaussian) Doppler broadening with HWHM Ai = 887. 



to decay as (using the original units) 



P 



-ft 



dAe- iA *p(A/ 7 )/ 7 



exp [--ft- A 2 t 2 /(41n2)] , 



(28) 



where we have used p(A/j) from Eq. (|25|) . Thus, losses 
induced by Gaussian broadening are non-Markovian. 
Since the time required for fast retrieval varies as i ~ 
1/(7^) ( see i f° r example, paper II), Doppler- induced 
losses become negligible compared to spontaneous emis- 
sion losses for sufficiently large d (d^> (Ai/7) 2 or, equiv- 
alently, d! ^ A1/7), and the system will behave as if 
there were no inhomogeneous broadening. The essential 
step in this last derivation is that the second moment 
(A 2 ) with respect to p(A/-f) is finite, and the result is 
thus applicable to any inhomogeneous profile falling off 
faster than a Lorentzian. In contrast, for a Lorentzian 
profile, the optical polarization decay would be 



P~exp(- 7 t- A^), 



(29) 



and we recover the effective rescaling of 7 up to 7 + Ai, 
which is equivalent to a simple rescaling of d down to d! . 

We now turn to the discussion of the optimal velocity- 
symmetric spin- wave modes S(z) obtained under Gaus- 
sian broadening using Eq. (|27p . These optimal modes 
are plotted (for forward retrieval) in Fig. [2] for 61 = 
0.17,0.67,3.69,14.25 and Ai = 887, which corresponds 
to the Rb Dl line at room temperature (assuming 27 is 
the natural linewidth). Indeed, at sufficiently high op- 
tical depth, the optimal mode approaches S(z) — y35 
(which we have derived in the homogeneously broadened 
case considered in paper II), since according to the ar- 
gument above, Doppler broadening plays no role at suf- 
ficiently high d. 

The above calculation (via Eq. (|27p) of the optimal 
spin wave yields the optimal retrieval efficiency. Using 
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the general time-reversal arguments presented in paper 
II, which still apply with Doppler broadening, we can, 
however, also use these optimal modes, to calculate the 
optimal efficiency for the combined process of storage 
followed by retrieval. The optimal symmetric modes 
(Sj(z) = S(z) for all j) for retrieval, which we have found 
above, are all real. Time reversal thus shows that the 
optimal storage into the symmetric mode, is obtained 
by time-reversing retrieval, and has the same efficiency 
as the retrieval. Note, however, that, in general, asym- 
metric modes may also be excited during storage. In 
this section, however, we assume that the atoms rether- 
malize during the storage time. This washes out any 
amplitude on asymmetric modes and the only relevant 
efficiency is the efficiency of storage onto the symmetric 
mode. The total maximum efficiency for storage followed 
by retrieval will thus be the square of the maximum re- 
trieval efficiency (obtained as the largest eigenvalue of the 
kernel in Eq. ([77]) ). With circles in Fig. [31 we show the 
maximum total efficiency for storage followed by back- 
ward retrieval for Aj = 887. The solid line in Fig. [3J is 
the Doppler- free maximum efficiency ^back(^) fr° m Fig- 4 
in paper II. The dotted line is the efficiency one would 
naively expect from a simple rescaling of the resonant op- 
tical depth, CackK), where d ' = d(7/Ai)V^E2. The 
dashed line is 5.8 (jr/^Ad 12 ) + 1/rf), which approximates 
the error fall-off of the points obtained from the full nu- 
merical optimization of Eq. ([27)) (circles) reasonably well 
and comes from the following heuristic model: t = l/("/d) 
is inserted into Eq. ([2"5)> . the exponential is expanded to 
first order, and a prefactor of 5.8 is introduced to match 
the Doppler-free error at large d, which is ~ 5.8 /d, as 
found in paper II. 

Above, we found the optimal spin-wave modes for re- 
trieval and the optimal retrieval efficiency by computing 
the eigenvector with the largest eigenvalue of the kernel 
in Eq. ([2T| . As in the homogeneously broadened case 
discussed in paper II, one can, in fact, interpret the it- 
erations used in finding the largest eigenvalue as the it- 
erations of retrieval followed by storage of time-reversed 
output and control pulses. To show this, we note that the 
retrieval in Eq. (fH))) inverse-Laplace-transformed back to 
t and the storage equation, which is obtained using the 
same steps, can be written as 

£out(t) = / d2m(t,z)S(l-z), (30) 
Jo 

S{z,f) = J dtff^T -t,z)e ia (£), (31) 

for some function to. Since these equations satisfy the 
general time-reversal form of Eqs. (35) and (36) in paper 
II, the same results as in the homogeneously broadened 
case apply. In particular, using Eqs. ([50]) and (pITj) , one 
can check that the maximization of retrieval efficiency 
through the iterative integration of the kernel in Eq. I|27p 
is equivalent to retrieval followed by time reversal and 
storage. From Eqs. ([50)) and ([3"T)) , it also follows that, 
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FIG. 3: Error 1 — 77 as a function of unbroadened optical depth 
d for different efficiencies r/. ??back(d) (solid) is the Doppler- 
free maximum total efficiency for any storage followed by 
backward retrieval. The maximum total efficiency for storage 
followed by backward retrieval with Doppler broadening of 
HWHM Ai = 887 (circles) does not follow <Tck(d') (dotted), 
which is the efficiency one would expect from naive rescaling 
of resonant optical depth, but follows more closely the heuris- 
tic model 5.8 (7r/(4d' 2 ) + 1/d) (dashed) described in the text. 



in order to optimize fast storage followed by fast forward 
retrieval, one should start with any trial input mode, 
store it, retrieve forward, time-reverse the whole process, 
and then iterate till convergence is reached, exactly as in 
the homogeneously broadened case of paper II. 

It is worth noting that a connection between time re- 
versal and optimal photon storage in the photon-echo 
technique was first made for the case of ideal, reversible 
storage in Refs. 0, dHHH. In the present paper and in 
Refs. [I], we extend this connection to a wide range 
of storage techniques in A-type media and show that op- 
timal storage is intimately connected with time reversal, 
even when the dynamics of the system are not completely 
reversible, and when the ideal unit efficiency cannot be 
achieved. 

We have thus shown that time-reversal iterations can 
be used to optimize storage followed by retrieval not only 
in the case of homogeneous broadened media discussed 
in paper I, paper II, and Ref. [l|, but also in the case of 
Doppler-broadened media. As explained in paper II, such 
time-reversal iterations not only constitute a convenient 
mathematical tool; they can, in fact, be used experimen- 
tally to find the optimal input modes. In particular, as 
explained in paper II, since the envelope £ of the quan- 
tum light mode obeys the same equations of motion as 
classical light pulses, one can first use the iterative pro- 
cedure to optimize the storage of classical light pulses, 
which can be easily measured and reversed, and then di- 
rectly apply this knowledge to the storage of quantum 
states of light. In fact, the first experiment on time- 
reversal-based optimization has already been successfully 
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carried out for classical light [13| and confirmed the va- 
lidity of the suggested procedure. 



III. INHOMOGENEOUS BROADENING 
WITHOUT REDISTRIBUTION BETWEEN 
FREQUENCY CLASSES DURING THE 
STORAGE TIME 

In the previous section, we treated the case when in- 
homogeneously broadened atoms redistribute themselves 
among different frequency classes during the storage 
time, which is the case in Doppler-broadened atomic va- 
pors for sufficiently long storage times. This redistribu- 
tion, however, does not take place in some other possi- 
ble experimental realizations, e.g., in Doppler-broadened 
atomic vapors with short storage times or in solid state 
media. Therefore, in this section, we consider what hap- 
pens when the redistribution among frequency classes 
does not take place. 

In the case of fast storage and retrieval, provided the 
7r pulse is applied at a sufficiently high power, it does not 
matter whether the two optical transitions are broadened 
independently or not, i.e., whether the \s) — \g) transition 
is broadened. However, in the case of adiabatic stor- 
age and retrieval the assumption that the control and 
the quantum field are always in two-photon resonance 
is crucial. Although the only regime of storage and re- 
trieval in inhomogencously broadened media we will con- 
sider in this section is the fast regime, the analysis will 
also be extendable to the adiabatic limit provided the 
\s) — \g) transition is homogeneously broadened. In this 
case, Eqs. ©-([7]) apply without modification. The proof 
in Appendix[A]that the retrieval efficiency is independent 
of the detuning and the control, therefore, also applies. 

Using the solution technique based on the Laplace 
transformation in time introduced in Sec. HH one can 
show that when storage is followed by forward retrieval 
and the inhomogeneous profile is Lorentzian, it actually 
does not matter whether the atoms redistribute them- 
selves among different frequency classes during the stor- 
age time or keep their frequencies unchanged: the same 
efficiency and output field are obtained. In this case, the 
results from paper II about homogeneous broadening are 
directly applicable if one replaces d by d! . This is, how- 
ever, not true for backward retrieval with a Lorentzian 
inhomogeneous profile or for retrieval in either direction 
with a Gaussian inhomogeneous profile. To obtain the 
efficiency in this situation, it is therefore necessary to 
take into account the fact that the transition frequency 
of each individual atom is the same during both storage 
and retrieval. Furthermore, by controlling and revers- 
ing the inhomogeneous broadening, one can even achieve 
rephasing of atomic excitations and, in fact, attain an in- 
crease in total efficiency relative to an unbroadened case 
7, 14]. An exhaustive study of the problem of storage 
followed by retrieval in media with no redistribution be- 
tween frequency classes during the storage time is beyond 



the scope of this paper. Here we restrict ourselves only to 
the investigation of fast storage followed by fast backward 
or forward retrieval in such media. We also include the 
possibility of reversing the inhomogeneous profile during 
the storage time as suggested in Refs. ^Ejj- I n particu- 
lar, in Sec. IIII Al we set up the equations for the problem 
of fast storage followed by fast retrieval in either direc- 
tion with and without the reversal of the inhomogeneous 
profile. In Sees. IIII Bl and IIII C[ we then discuss the re- 
sults that these equations yield for the cases of storage 
followed by backward retrieval without and with the re- 
versal of the inhomogeneous profile, respectively. 

A. Setup and solution 

In this section, assuming that the redistribution be- 
tween frequency classes takes place neither during the 
processes of storage and retrieval nor during the storage 
time, we set up and solve the problem of fast storage 
followed by fast retrieval in the forward or backward di- 
rection with or without the reversal of the inhomogeneous 
broadening. Any storage with no redistribution between 
frequency classes (not only the fast limit) can be com- 
puted numerically with discrete frequency classes for any 
kind of inhomogeneous profile and any control. To do 
this, one can just use Eqs. ©-([7]) for both storage and 
retrieval and, depending on the direction of retrieval and 
on whether the inhomogeneous profile is reversed during 
the storage time, make the appropriate modification to 
the stored spin waves Sj(z,T) before retrieving it. How- 
ever, as we now show, the case of fast storage and fast 
retrieval can also be solved almost completely analyti- 
cally using Laplace transformation in time, t —> v [15j |. 
Before the storage tt pulse is applied, f2 = in Eq. ©. 
Then for the fast storage of a resonant input mode £- m (t), 
Eqs. ([5]) and I© can be solved to give 

P {~z,v) = iVd ^ £ in (v)e- dS ^ v \ (32) 

1 + v + iAj 

where v is the Laplace variable, and f(v) is defined in 
Eq. (|20p and is computed for various inhomogeneous pro- 
files in Eqs. (|2T|) - ([2H1) . 

To find the initial conditions for the subsequent re- 
trieval, we take the inverse Laplace transform u — > t = T 
of Eq. (f3"2"]l and multiply Pj by —1 to account for the 
two 7r pulses (i.e., the storage and retrieval tt pulses). If 
we are interested in backward retrieval, Pj(z) should be 
flipped to Pj(l — z). If we are interested in reversing 
the inhomogeneous profile, the frequency classes should 
be reversed. Thus, for example, for backward retrieval 
with the reversal of inhomogeneous broadening, the ini- 
tial condition for retrieval is Pj(z, T r ) = —P-j(l — z, T). 
Using Eqs. ([5]) and ^ with f2 = 0, we can then imple- 
ment fast retrieval. The time Laplace transform of the 
output field can then be found to be equal to 

£out{v) = £- 1 {A{v 1 v')B(v, v')E- m {v')} v ,^ f , (33) 
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where Cr 1 indicates that we should take the inverse 
Laplace transform v' — ► T, A(v, v') depends on the di- 
rection of retrieval and is given by 



A(v,v') 
A(v,v') 



e -d(f(v)+f(v')) _ 1 

f(v) + f(v r ) ' 

e ~df(v) _ e -df(v') 

f(v) - W) 



(34) 
(35) 



for backward and forward retrieval, respectively, and 
B(v,v') depends on whether the inhomogeneous profile 
is reversed or not and is given by 



B(v,v') = 
B(v,v') = 



M + f(y') 

2 + v + v' ' 



(36) 
(37) 



for reversed and not reversed cases, respectively. 

For homogeneous broadening (Eqs. (|21[) and (HU)) and 
for a Lorentzian inhomogeneous profile (Eqs. (|23|) and 
(I24p ). the inverse Laplace transforms v' — > T and v — > t 
can be taken analytically in terms of Bessel functions 
and convolutions. The case of homogeneous broadening 
has been studied in paper II, while the analytical answer 
for the Lorentzian case is too complicated to yield any 
significant insight. For a Gaussian inhomogeneous profile 
(Eqs. (|23|) and (f!?rj|) ), the inverse Laplace transforms have 
to be taken numerically. In all three cases, the efficiency 
can be computed via Eq. (fl"4)) without taking the v — > t 
inverse Laplace transform. 

If we are interested not only in computing the total 
efficiency of storage followed by retrieval for some given 
input photon mode, but also in maximizing the efficiency 
with respect to the input mode shape, we can again take 
advantage of time reversal. We will show now that in 
all four cases (i.e., either of the two retrieval directions, 
with or without the reversal of broadening) , the optimal 
input pulse shape can be found by starting with any trial 
input, carrying out storage followed by retrieval, then 
time-reversing the output, and iterating the procedure 
till convergence is reached. To begin the proof, we note 
that A(v,v') and B(v,v') are symmetric with respect to 
the exchange of the two arguments. Therefore, Eq. (f33|) 
and the convolution theorem for Laplace transforms im- 
ply that we can write 



£out(£) = / dt'£ in (t')m'(t, f - t') 



(38) 



for some function m' that is symmetric with respect to 
the exchange of its two arguments. One can also check 
that m', and hence the optimal input mode, are real. 
Assuming, therefore, a real £i n {i), the total efficiency is 



where the kernel k to t(t,t') (the subscript "tot" stands for 
the total efficiency, i.e., storage followed by retrieval) is 
given by 

fetat (*,?)=/ dt"m'(t",f-t)m'(t",f-t / ), (40) 
Jo 

where we assumed T is sufficiently large that the interval 
[0,T] includes the whole retrieved pulse. To find the 
optimal £- m (t), one can thus start with a trial input mode 
£\ (t) and iterate the action of the kernel according to 



di£i(t)fctot(*,?). 



(41) 



Using the symmetry of m! , one can immediately see that 
this iteration is equivalent to carrying out storage fol- 
lowed by retrieval, time-reversing the output, repeating 
the procedure, and time-reversing the output again. We 
will use these time-reversal iterations in the following sec- 
tions to compute the optimal input modes. It is impor- 
tant to note that, as we have explained at the end of 
Sec. IIIB1 the time-reversal iterations that we have just 
described not only constitute a convenient mathematical 
tool; they can, in fact, be used experimentally to find the 
optimal input modes. 

In this section, we set up and solved the problem of fast 
storage followed by fast backward or forward retrieval 
with or without the reversal of the inhomogeneous profile. 
In the next two sections, we would like to analyze these 
solutions for the case of backward retrieval. 



Storage followed by backward retrieval 



7J = 



(It 



dt'& m (t)£ in (t')k tot (t, t') , (39) 



In this section, using the results of Sec. lIII Al we would 
like to study the efficiency of fast storage followed by fast 
backward retrieval in inhomogeneously broadened me- 
dia without redistribution between velocity classes, and 
without the reversal of the inhomogeneous broadening 
during the storage time. This problem is motivated by 
solid state implementations such as the one described in 
Refs. [H, [l6| , where the line shape is created by pump- 
ing back some atoms from a broad spectral hole feature. 
In this situation, one can consider what happens when 
one expands the spectral region and pumps more and 
more atoms into the absorptive feature. In this case, 
the resonant optical depth d' is expected to be indepen- 
dent of the width of the spectral feature that is pumped 
back into the absorption profile, whereas the unbroad- 
ened optical depth d would increase with an increasing 
number of atoms being pumped back. From the analysis 
of Sec. UU we expect the behavior to be different depend- 
ing on whether the line shape falls off as a Lorentzian 
or faster. In particular, we found in Sec. [Til that for a 
well-localized inhomogeneous line (e.g., a Gaussian), the 
losses due to dephasing are non-Markovian and, there- 
fore, scale as 1/d' 2 , in contrast to the losses due to the 
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exponential polarization decay with the rate 7, which 
scale as 1/d. As another application of this idea, we show 
in this section that, for a Gaussian inhomogeneous pro- 
file, by increasing d, the error during fast storage followed 
by fast backward retrieval can be lowered from the 1 j d! 
scaling for the homogeneously broadened line (d = d') 
to the 1/d' 2 scaling for sufficiently large d. Experimen- 
tally, d! can be quite large (d! ss 50 should be feasible 
by using a sufficiently high impurity concentration and a 
sufficiently long sample length so that this could 

potentially yield a very high efficiency. 

We begin the analysis by using the time reversal iter- 
ations suggested at the end of Sec. IIII Al to compute the 
optimal input pulse and the maximum possible efficiency 
at various values of d and d! assuming the inhomogeneous 
profile is Gaussian (Eq. ([25)) ). To compute the standard 
deviation a of the Gaussian profile from the values of d 
and d! , we use the following relationship: 



d! = dx 



dAp(A) 



1 



1 + A 2 



e 2 " 

2 a 



"erfc 



Via 



where the second factor in the integrand is the homoge- 
neous line shape of HWHM equal to 1 (recall that our 
frequencies are rescaled by 7). 

In Fig.|4j we show (solid lines) the optimal input modes 
for d' = 20 and d = 20, 60, and 120. We see that these 
optimal input modes (for this value of d') are very sim- 
ilar to the optimal input modes (dashed lines) for the 
same d but without the inhomogeneous broadening (i.e., 
d' = d). The optimal modes thus have a duration of ap- 
proximately l/(dy), as discussed in Sec. VII of paper II, 
and consist of a roughly triangular pulse preceded by a 
few "wiggles." These wiggles can be traced back to the 
zeroth-order Bessel function of the first kind (i.e., Jo) in 
Eq. (48) in paper II and can be associated with Rabi 
oscillations between the electric field £ and the optical 
polarization P. In Fig. |4j we assume that the storage 
7r pulse is applied at t = 0. Although the true optimal 
input pulses extend to t = —00, they can, in practice, 
be truncated after about two or three wiggles without 
decreasing the efficiency by more than 10~ 4 . 

To verify the prediction that, for sufficiently large d, 
the error should be limited by 1/d' 2 , we computed the 
optimal (smallest) error by optimizing with respect to 
the input mode at different values of d' and d. The op- 
timization was done numerically using the time-reversal 
iterations suggested at the end of the previous section. 
Figure 03a) shows a log-log surface plot of this optimal 
(smallest) error 1 — 77 as a function of d' and d/d! . As 
expected, we find that, for any fixed d! , the error is very 
well approximated by c\ +C2/ 1 d at large d, where the con- 
stants ci(d') and C2(d') depend on d' . C2(d') is of order 
unity and increases approximately linearly from about 
0.2 at d' = 2 to 4.2 at d! = 20. a(d') represents the 
(i'-limited error, i.e., the limit d/d' — > 00, when the 1/d 
error becomes negligible. This d'-limited error can be 
seen at the d/d' = 10 3 edge of the box in Fig. EDJa) and is 




FIG. 4: The solid lines show the optimal input modes for d' = 
20 and d = 20, 60, and 120 (the value of d is indicated on the 
plot) for storage followed by backward retrieval in a medium 
with a Gaussian inhomogeneous profile, without the reversal 
of the inhomogeneous broadening, and without averaging over 
frequency classes during the storage time. The dashed lines 
show the corresponding optimal modes for d' — d, i.e., for the 
case of no inhomogeneous broadening. 



also plotted separately as a thin solid line in the log-log 
plot in Fig.OJb). The dotted line in Figs.[5fa) and[5]Jb) is 
14.2/ d' 2 and is shown to indicate that the d'-limited error 
indeed scales as 1/d' 2 . To include the requirement that 
the efficiency drops to zero at d! = and to reproduce 
the 14.2/d' 2 dependence that we see at larger d' , a dashed 
heuristic curve l/(l + 14.2/d' 2 ) is shown in Fig. |5fb) . We 
see that it matches the exact value of the e?'-limited error 
(thin solid line) very well. 

We compare this e?'-limited error to the smallest pos- 
sible error for the homogeneously broadened case d = d! , 
which can be seen at the d/d' = 1 edge of the box in 
Fig. 02a) and which is also plotted separately as a thick 
solid line in Fig. [5^b). This thick solid line is the same 
as the solid curve in Fig. 4 of paper II and in Fig. [3] of 
the present paper. Since we know from the discussion 
in paper II that at high enough values of d (= d'), this 
error scales as ~ 5.8/d, we plot this 5.8/d scaling as the 
dash-dotted line in Figs. (Ha) andOJb). We thus see that, 
for a fixed d! , by increasing the number of atoms in the 
absorption line so that we go from d/d' = 1 to d/d' — ► 00 
(in such a way that the inhomogeneous line is Gaussian) , 
one can significantly lower the optimal (smallest) error 
from a 1/d' homogeneous error to a 1/d' 2 inhomogeneous- 
broadening-limited error. Although we have investigated 
only backward retrieval, we expect that for the case of 
optimal storage followed by forward retrieval, the same 
result will apply, and that the optimal 1 / d! homogeneous 
error can also be reduced to a d'-limited 1/d' 2 error. We 
also expect these error scalings to hold not only for the 
Gaussian inhomogeneous profile but also for any inhomo- 
geneous profile (such as, for example, a square profile) 
whose tails fall off faster than Lorentzian. 
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FIG. 5: (Color online) (a) As a function of d! and d/d! , the 
optimal (smallest) error of fast storage followed by fast back- 
ward retrieval (surface). 5.8/d (dash-dotted line) is the lim- 
iting behavior, found in paper II, of the homogeneous error 
(d/d' = 1) for large enough d (= d!). 14.2/ d' 2 (dotted) ap- 
proximates the d'-limited error (d/d 1 — > oo) near d! ~ 10 — 20 
and, probably, above, (b) As a function of d! , dash-dotted 
and dotted lines are the same as in (a). The thick and thin 
solid lines are the d/d' — 1 and d/d' — > oo error curves, re- 
spectively, which can be read out from the surface plot in (a). 
1/(1 + 14. 2/d' 2 ) (dashed line) is a heuristic curve that matches 
the d'-limited error (thin solid line) very well. 



To summarize, we have shown that, even without any 
additional experimental requirements, the storage and re- 
trieval by two fast 7r pulses may be a very promising route 
to a quantum memory. Indeed, if the absorptive fea- 
ture is sufficiently well localized (falling off faster than a 
Lorentzian), the error may be limited to ~ 1/d' 2 . In a 
practical realization, it would, however, be necessary to 
include other imperfections such as imperfect n pulses, 



imperfect synchronization between the storage ir pulse 
and the input pulse, and the limitations associated with 
the creation of an absorption line. Furthermore, the op- 
timal scenario described in this section applies only to a 
single input pulse shape at any given values of d and d' . 
For other pulse shapes, it may be advantageous to com- 
bine the ideas of this section with the reversible broad- 
ening of the next section and optimize with respect to 
both the width (and shape) of the original nonreversible 
line and the width (and shape) of the extra reversible 
broadening. 

C. Storage followed by backward retrieval with the 
reversal of the inhomogeneous profile 

In the previous section, we assumed that the inhomo- 
geneous broadening had a fixed distribution, e.g., due 
to different environments for each individual atom. In 
this section, using the results of Sec. lIII Al we investigate 
the possibility of improving the efficiency of fast stor- 
age followed by fast backward retrieval by adding and 
reversing inhomogeneous broadening in an originally ho- 
mogeneously broadened medium. 

The possible advantage of inhomogeneous broadening 
was first considered in Ref. where it was noted that 
Doppler broadening automatically reverses during back- 
ward retrieval. This results in a reversal of the dephasing 
occurring during storage and gives rise to photon echo. 
In Refs. @, [l4| it was then suggested that similar ef- 
fects could be realized in solid state systems. Under the 
name of controlled reversible inhomogeneous broadening 
(CRIB), the authors of Ref. Q suggest implementing the 
equivalent of the fast storage protocol considered in pa- 
per II, but in addition they suggest controllably adding 
inhomogeneous broadening to the transition and then re- 
verse the broadening during retrieval to obtain a rephas- 
ing [ItJ ■ Several experimental groups are currently work- 
ing on the realization of CRIB [EMEU- We will show 
that, although the introduction of reversible inhomoge- 
neous broadening can improve the efficiency of fast stor- 
age of a single pulse, the improvement relative to the 
fast storage technique without inhomogeneous broaden- 
ing is limited. We will also show that CRIB can perform 
slightly better than optimal adiabatic storage in a ho- 
mogeneously broadened medium (discussed in Sees. VI 
B and VI C of paper II), but only for short pulses, for 
which the adiabaticity condition Td*y ^ 1 is not satisfied. 

For concreteness, we consider the storage of a resonant 
Gaussian-like pulse of variable time duration T defined, 
as in Eq. (40) in paper II, by 

S in (t) = A{e -mi/f-o.tf _ e -T.sy^f (43) 

and shown in Fig. 3 of paper II. A « 2.09 is a normaliza- 
tion constant. We also restrict ourselves to the situation 
where the medium is initially homogeneously broadened. 
Although this is often a good approximation, the hole- 
burning technique in the solid state [16] , for example, will 
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always result in some residual inhomogeneous broaden- 
ing of the prepared line, something that one may have to 
take into account in a complete assessment of the perfor- 
mance of CRIB. 

To investigate the performance of CRIB, we use the 
techniques described in Sec. IIII Al Let us initially as- 
sume Tj < 1, so that the decay can be ignored. As 
a test case, we take the resonant Gaussian-like pulse of 
Eq. and implement fast storage of it (with a tt pulse 
at t = T) followed by fast backward retrieval. From the 
discussion of CRIB [7], it is expected that adding some 
broadening and thereby increasing the width of the ab- 
sorption line, at the expense of a decreased resonant op- 
tical depth, may increase the total efficiency Excessive 
broadening will, however, make the medium transparent 
and decrease the efficiency. At each value of Tdj, there 
is thus an optimal inhomogeneous width. In Fig. [51 the 
total efficiency is plotted as a function of Td'y, for the ho- 
mogeneously broadened case (dash-dotted line) and for 
an inhomogeneously broadened medium (solid lines) with 
Gaussian (G) or Lorentzian (L) inhomogeneous profiles 
optimized with respect to the inhomogeneous width. The 
horizontal dashed line is 1. These curves are calculated 
by numerically computing, at each Td'y, the efficiency as 
a function of the inhomogeneous width and then finding 
the width that gives the maximum efficiency at the cho- 
sen value of Tdj. Note that, even though we neglect the 
decay (7T < 1), the quantity Td-y = g 2 NTL/c may still 
attain a non-negligible value for a sufficiently high optical 
depth d> 1. In other words, the curves in Fig. [6] repre- 
sent the limit T7 — ► and d — ► 00 with finite Td'y. Below 
Td'y ~ 1, adding broadening only lowers the efficiency, 
so that the optimal curves (solid) join the unbroadened 
curve (dash-dotted). We see, however, that at intermedi- 
ate values of Td'y ~ 10, introducing reversible inhomoge- 
neous broadening can increase the efficiency. The gain in 
efficiency is, however, limited. In an experimental real- 
ization, one should therefore evaluate whether this gain 
justifies the additional experimental efforts needed to im- 
plement the controlled reversible inhomogeneous broad- 
ening. 

In Fig. [71 we plot A\T as a fuction of Td'y, where Ai is 
the optimal HWHM of the Gaussian (G) or Lorentzian 
(L) inhomogeneous profile used to construct Fig. The 
points are connected with straight lines for better visi- 
bility. The dashed and dotted lines are 1.4(Td 7 - 2) 1 / 2 
and 2.25(Te?7 — 2) 1 / 2 , respectively, which indicates that 
the optimal inhomogeneous HWHM Ai scales approx- 
imately as y^Tdy/T, which is different from the naive 
guess Ai ~ 1/T. We also see that, at a given Td'y, the 
optimal Gaussian profile is wider; and, moreover, from 
Fig. [SJ we see that optimal Gaussian broadening gives 
greater efficiency than optimal Lorentzian broadening. 
One could imagine that these two results are the con- 
sequence of the Gaussian frequency profile of the input 
pulse. However, we ran an equivalent simulation with 
an input pulse that has a Lorentzian spectrum; and also 
in this case the optimal Gaussian profile is wider than 
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FIG. 6: Total efficiency of fast storage followed by fast back- 
ward retrieval as a function of Td'y in the limit when Ty — > 0, 
d — > 00, but Tdy is finite. The curves show the efficiency with- 
out inhomogeneous broadening (dash-dotted line) and with 
(solid lines) Gaussian (G) or Lorentzian (L) inhomogeneous 
profiles optimized with respect to the inhomogeneous width. 
The horizontal dashed line is 1. 



the optimal Lorentzian profile and the optimal Gaussian 
broadening gives a greater efficiency than the optimal 
Lorentzian broadening. These results reflect the fact that 
the storage we are considering is a dynamical process and 
is therefore not accurately described by its continuous- 
wave absorption profile. We believe that the advantage 
of Gaussian broadening over Lorentzian comes from the 
fact that, as shown in Sec.[TTl due to their non-Markovian 
nature, the dephasing losses associated with Gaussian 
(and, hence, fast-falling) inhomogeneous broadening are 
smaller than the losses associated with Lorentzian broad- 
ening. 

To investigate the performance of the storage and re- 
trieval protocol at a finite optical depth, we now relax 
the assumption 7T <C 1. With the pulse considered here 




Tdy 



FIG. 7: AiT as a function of Tdy, where Ai is the optimal 
HWHM of the Gaussian (G) or Lorentzian (L) inhomogeneous 
profile. The points are connected with straight lines to guide 
the eye. The dashed and dotted lines are 1.4(Td 7 -2) 1/2 and 
2.25(Td 7 - 2) 1/2 , respectively. 
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(Eq. (|43f ). the effect of the spontaneous emission 7 on 
fast storage with CRIB can be estimated based on simple 
arguments: the pulse is symmetric around T/2 so the ex- 
citation spends on average a time T/2 in the sample both 
during storage and retrieval. In each of these processes, 
the efficiency is therefore decreased by [exp(— 7T/2)] 2 
so that the total efficiency is reduced by approximately 
exp(— 27T). We have explicitly verified this simple esti- 
mate for a few cases and found it to be true both with 
and without broadening. The optimization of broadening 
without decay therefore gives the same optimal inhomo- 
geneous width as with decay. Since we have in Sec. VI D 
of paper II calculated adiabatic efficiencies for the same 
pulse shape, we can now compare the performance of fast 
storage with and without CRIB to adiabatic storage. In 
Fig. [5Ja), with d = 100, we compare the storage of the 
Gaussian-like pulse of duration T of Eq. (|43|) (shown in 
Fig. 3 of paper II) followed by backward retrieval using 
"optimal" adiabatic storage (dotted) or using fast stor- 
age without inhomogeneous broadening (dash dotted) or 
with reversible optimal- width Lorentzian (L) or Gaussian 
(G) broadening (solid). The horizontal dashed line is the 
optimal adiabatic efficiency, while the second dashed line 
is exp(— 27T), by which fast efficiencies are rescaled rel- 
ative to the T7 — > limit of Fig. O Figure [2Kb) is the 
same as Fig.(5Ja) but for d = 1000. 

First of all, from Fig. [5J it is clear that, when Tj > 1, 
fast storage efficiency, with or without CRIB, deterio- 
rates because the spontaneous emission decreases the to- 
tal efficiency by exp(— 27T). Optimal adiabatic storage, 
on the other hand, does well in this limit (provided we 
have a reasonable d). Moreover, in the adiabatic limit 
Tdj 3> 1, optimal adiabatic storage with homogeneous 
broadening is always more efficient than fast storage with 
or without CRIB. This follows from the fact that in 
the adiabatic limit (Tdj 3> 1) the error in fast storage 
(w 1 — cxp(— 27T) sa 27T) is greater than the error in 
adiabatic storage (~ 5.8/d). 

Secondly, when Tdj < 1, neither of the methods does 
very well. This follows directly from the general time- 
reversal argument described in detail in paper II. Ac- 
cording to these arguments, the optimal storage is ob- 
tained as the time reverse of optimal retrieval. One can- 
not therefore store faster than one can retrieve. The 
fastest retrieval is obtained by using the fast retrieval 
method where the excitation is transferred from state |s) 
into state |e) with a n pulse. With this procedure, all 
atoms radiate in phase and, by constructive interference, 
give a short output pulse of duration T ~ l/{drf), as 
explained in Sec. VII of paper II. This procedure gives 
the fastest possible retrieval, and its time reverse is the 
fastest possible storage (which works optimally only for 
certain input mode shapes). The storage and retrieval 
of any mode thus becomes inefficient for Td^ < 1. The 
particular Gaussian mode function that we consider here 
does not correspond to the optimal mode and therefore 
its fast storage-plus-retrieval efficiency (dash-dotted line 
in Fig. [5]) does not reach the optimal efficiency (horizon- 
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FIG. 8: Storage and backward retrieval with and without 
reversible inhomogeneous broadening at finite optical depth 
d = 100 (a), and (b) 1000. The curves show the total efficiency 
of storage followed by backward retrieval of the resonant 
Gaussian-like pulse of duration T of Eq. (143 [) (shown in Fig. 3 
of paper II) with "optimal" resonant adiabatic storage (dot- 
ted) or with fast storage without inhomogeneous broadening 
(dash dotted) or with reversible optimal-width Lorentzian (L) 
or Gaussian (G) broadening (solid). The horizontal dashed 
line is the optimal adiabatic efficiency while the second dashed 
line is exp(— 27T), by which fast efficiencies are rescaled rel- 
ative to the T7 -> limit of Fig. |6] 



tal dashed line in Fig. [5J. 

Third, when Tdj < 25, reversible inhomogeneous 
broadening (CRIB) does help, and, with it, the fast 
method may do slightly better than the adiabatic method 
without inhomogeneous broadening. An interesting pos- 
sibility is whether the controlled addition and reversing of 
inhomogeneous broadening could improve the adiabatic 
storage in the regime when the adiabaticity condition 
Td'y 3> 1 is not satisfied, but this investigation is beyond 
the scope of the present paper. 

To summarize this investigation of the possible advan- 
tages of introducing a reversible inhomogeneous broad- 
ening, we conclude that it does provide an improvement 
of the efficiency, but that this improvement is limited. 
The fast storage technique where one just applies a res- 
onant 7T pulse without the additional broadening gives 
comparable results for the pulse shape that we have con- 
sidered here. Moreover, even with an optimized inho- 
mogeneous width, fast storage with CRIB performs only 
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slightly better than optimal adiabatic storage in the origi- 
nal homogeneously broadened system and only when the 
input pulse does not satisfy the adiabaticity condition 
Tdj 3> 1; while, for pulses satisfying the adiabaticity 
condition, adiabatic storage performs better. Intuitively 
one might expect that a homogeneously broadened ab- 
sorption line would not be able to efficiently store an 
input pulse with a bandwidth that is much longer than 
the width of this line. One could have therefore expected 
a large efficiency gain from the use of CRIB to shape 
the atomic line to match the spectral profile of the input 
photon wave packet. The reason why such line shaping 
is not necessarily much more effective than using the un- 
modified line is because storage is a dynamical process; 
therefore, the relationship between the storage capabil- 
ity of the medium and its continuous wave absorption 
spectrum is not trivial. 

We should add, however, that we have here only con- 
sidered a specific input shape and the picture may be 
different if other inputs are considered. We also empha- 
size that, unlike most results in this paper and in papers 
I and II, the efficiencies presented in this section do not 
represent the true optimum for our input pulse, since 
one could optimize, for example, the time, at which the 
storage tt pulse is applied. Finally, it is worth noting 
that if one is free to choose the shape and duration of 
the input pulse, then for any given d and any given in- 
homongeous profile, as explained in Sec. IIII Al one can 
use the time-reversal iterations to find the optimal input 
pulse and the maximum efficiency. We have checked for 
a Gaussian inhomogeneous profile that at any given d, as 
the inhomogeneous width increases from zero, the maxi- 
mum efficiency drops. Although we have not checked this 
statement for other inhomogeneous profiles, we believe it 
to be generally true that, if one has the freedom of op- 
timizing the shape and duration of the input pulse, then 
the addition and subsequent reversal of inhomogeneous 
broadening only lowers the maximum efficiency. 

We note, however, that our result that the adiabatic 
storage efficiency is always comaparable to or better than 
the efficiency of fast storage (with or without CRIB) 
has been found assuming that in the adiabatic storage 
the \g) — \s) transition is not inhomogeneously broad- 
ened. Although this holds in Doppler-broadened atomic 
vapors, the broadening of the \g) — \s) transition may 
be hard to suppress in many of the systems considered 
for photon storage with CRIB (such as the rare-earth- 
ion-doped crystal [lj]). Therefore, although adiabatic 
storage in some of these systems seems possible [l8j], 
fast storage may still be a better option. Moreover, al- 
though we showed that, for fast storage of a single input 
mode, adding inhomogeneous broadening may provide 
only small gains in efficiency, these gains might be much 
more significant if several time-separated modes are to 
be stored together and the shape of the inhomogeneous 
profile is allowed to be optimized. Further investigations 
are required to clarify these issues. 



IV. SUMMARY 

In conclusion, we have extended in this paper the anal- 
ysis of photon storage in papers I and II to include the 
effects of inhomogeneous broadening. In particular, we 
showed that in Doppler-broadened atomic vapors, at high 
enough optical depth, all atoms contribute coherently as 
if the medium were homogeneously broadened. We also 
showed that high-efficiency photon storage (error scaling 
as ~ 1/d' 2 , where d! is the observed optical depth) can 
be achieved in solid state systems by creating a station- 
ary spectrally well-localized inhomogeneous profile. Fi- 
nally, we demonstrated that the addition of reversible in- 
homogeneous broadening (CRIB) to an originally homo- 
geneously broadened line does provide an improvement 
of the efficiency of fast storage followed by fast retrieval, 
but that this improvement is limited: in particular, in 
the adiabatic limit Tdj ^> 1 optimal adiabatic storage 
outperforms fast storage with or without CRIB. These 
results aim at understanding the fundamtental limits for 
storage imposed by the optical depth of the medium. For 
a complete investigation of the photon storage problem 
in inhomogeneously broadened media, there are several 
other effects and experimental imperfections that should 
be included, such as, for example, velocity changing col- 
lisions during the processes of storage and retrieval in 
Doppler-broadened gases and imperfect synchronization 
between the input pulse and the storage control pulse. A 
study of the former is in progress. 

The presented optimization of the storage and retrieval 
processes in inhomogeneously broadened media leads to a 
substantial increase in the memory efficiency. We there- 
fore expect this work to be important in improving the ef- 
ficiencies in current experiments, where the optical depth 
is limited and where inhomogeneous broadening plays an 
important role, such as in Doppler-broadened atoms in 
warm vapor cells [Jland in inhomogeneously broadened 
solid state samples p, [l4| . 

Note added. Recently, a related paper appeared [li|, 
which discusses some of the issues considered in this pa- 
per. 
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APPENDIX A: INDEPENDENCE OF 
FREE-SPACE RETRIEVAL EFFICIENCY FROM 
CONTROL AND DETUNING 



where L 1 stands for inverse Laplace transforms u 
and u' — > z! with the evaluation of both at z — z' 
From Eqs. (|A2j) and JA3|, it follows that 



We showed in Sec. Ill Bl that for the case of inhomo- 
geneous broadening that leaves the \s) — \g) transition 
homogeneously broadened, the retrieval efficiency in the 
cavity model is independent of the control field and the 
detuning (provided all excitations are pumped out of the 
atoms). In this appendix, we show the the same result 
holds in the free-space model. 

We consider a general situation in which Sj(z,t = 0) 
are not necessarily equal for different velocity classes j 
(we assume here that retrieval begins at t = 0). Laplace 
transforming Eqs. ([5])-([7]) in space (z — ► u) and using 
£\ n = 0, we obtain 



£ = 



.Vd 
i — P, 
u 



(Al) 



dtPj 
dtSj 



-(1 + z(A + Aj))Pj - ~^p]P + i nS 3 ,(A2) 



(A3) 



where P = ^ fc ^/pkP k - The retrieval efficiency is then 
— / dtP 3 {u,t) [P k (u'*,t) 

UU J Q 



(A4) 



i (Pjiu, t) [P k (u'*,t)]*+ S^u, t) [S k (u'*,t)] ' 

-(2 + *(Aj - A k ))P,(u,t) [P k (u'*,t)]* (A5) 

i^pjp(u,t) [p^^Y-^p^t) [p(u'*,t)y. 



If M is the number of velocity classes, Eq. (|A5[) stands for 
M 2 equations in M 2 variables P%P k - We can write them 
in a matrix form and, in principle, invert the M 2 x M 2 
matrix on the right-hand side and thus solve for PjP k as 
a linear combination of 4;(P a P£ + S a S^) for various a 
and b. Inserting this into Eq. (|A4|) . applying the funda- 
mental theorem of calculus, and assuming the retrieval is 
complete (i.e., no excitations remain in the atoms), the 
retrieval efficiency can be expressed as C^ 1 of a linear 
combination (with u- and it'-dependent coefficients) of 
S^i^O) [Sb(u'* , 0)]*, and is thus independent of control 
and detuning. 
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